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Abstract 

Purpose: To design a mathematical model that can predict the relationship between the ganglion cell complex (GCC) 
thicl<ness and visual field sensitivity (VFS) in glaucoma patients. 

Design: Retrospective cross-sectional case series. 

Method: W\Xhm 3 months from VFS measurements by the Humphrey field analyzer 10-2 program, 83 eyes underwent 
macular GCC thickness measurements by spectral-domain optical coherence tomography (SD-OCT). Data were used to 
construct a multiple logistic model that depicted the relationship between the explanatory variables (GCC thickness, age, 
sex, and spherical equivalent of refractive errors) determined by a regression analysis and the mean VFS corresponding to 
the SD-OCT scanned area. Analyses were performed in half or 8 segmented local areas as well as in whole scanned areas. A 
simple logistic model that included GCC thickness as the single explanatory variable was also constructed. The ability of the 
logistic models to depict the real GCC thickness/VFS in SAP distribution was analyzed by the test of goodness-of-fit. The 
significance of the model effect was analyzed by analysis of variance (ANOVA). 

Results: ScatXex plots between the GCC thickness and the mean VFS showed sigmoid curves. The test of goodness-of-fit 
revealed that the multiple logistic models showed a good fit for the real GCC thickness/VFS distribution in all areas except 
the nasal-inferior-outer area. ANOVA revealed that all of the multiple logistic models significantly predicted the VFS based 
on the explanatory variables. Although simple logistic models also exhibited significant VFS predictability based on the GCC 
thickness, the model effect was less than that observed for the multiple logistic models. 

Conclusions: The currently proposed logistic models are useful methods for depicting relationships between the 
explanatory variables, including the GCC thickness, and the mean VFS in glaucoma patients. 
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introduction 

Glaucoma is a neurodegenerative disease associated with the 
progressive loss of retinal ganglion cells (RGCs). Visual field testing 
has shown that loss of RGCs leads to decreased sensitivity to light 
stimulation. The current structure-function relationship hypothesis 
states that there is a retinal ganglion cell (RGC) functional reserve 
or RGC redundancy in early clinical glaucoma. This hypothesis is 
consistent with the histological finding that as many as half of the 
RGCs and their axons may be lost before there is detectable visual 
function loss, as measured by standard automated perimetry (SAP) 



[1,2]. Therefore, the loss of RGCs likely leads to the atrophy of the 
ganglion cell layer (GCL). 

Newer versions of optical coherence tomography (OCT) that 
incorporate spectral-domain (SD) technology provide higher scan 
resolutions and scan speeds. It has been shown that this technique 
can measure the thickness of the three or four innermost layers in 
the macula that are atrophic in glaucoma [3]. A further OCT 
study additionally confirmed there was a reduced thickness in the 
innermost three to four layers in the retina of glaucoma patients 
[4]. Upon development of SD-OCT, one of the first areas 
measured was the thickness of the ganglion cell complex (GCC), 
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which consists of the retinal nerve fiber layer, GCL, and the inner 
plexiform layer (IPL) [5]. 

The continued development and evolution of the structure 
measurement methods have assisted researchers in their attempts 
to link the structural measurements to the visual functions. 
Previous studies that examined the structure-function relationship 
in glaucoma have developed one-to-one models between the 
structure (e.g., GCC thickness) and function (visual field sensitivity) 
[6-11]. However, other parameters such as age or spherical 
equivalent may also have a significant relationship with the 
structure measurements, such as the NFL or GCC thicknesses. In 
addition, many of these previous studies were not able to report 
definitive findings due to limited patient backgrounds, and because 
these studies only evaluated a small number of parameters. 
Therefore, to expand and clarify our understanding of the 
structure-function relationship in glaucoma, multiple regression 
analysis of previously excluded parameters will need to be 
performed. The purpose of current study was to design a 
mathematical model that cotdd use the GCC thickness and other 
parameters, such as age, sex, or spherical equivalent, to predict the 
SAP-determined visual field sensitivity (VFS). 

The relationship between the SAP-determined VFS and the 
GCC thickness does not appear to be a simple linear pattern. VFS 
reaches a plateau and does not increase without limitation, even 
when the GCC is sufficiently thick. On the other hand, steep 
decreases in the VFS seem to occur after the GCC thickness falls 
below a certain threshold. In addition, regardless of the type of 
SAP used, there are minimum and maximum detectable limits of 
VFS. Because of these characteristics and restrictions, plots of the 
GCC thickness/VFS in SAP curves result in so-called sigmoid or 
z-shaped curves. Thus, these curves cannot be represented by a 
polynomial formtjla. Furthermore, with the exception of the GCC 
thickness, the effects of the other variants of VFS have yet to be 
definitively clarified. Therefore, the use of a nonlinear logistic 
analysis may be appropriate when examining such a complex 
pattern. 

In this study, we analyzed the relationship between the local 
GCC thickness measured by SD-OCT and the VFS measur(^d by 
SAP in eyes with glaucoma in conjunction with other explanatory 
variables determined by logistic regression analyses. We also 
evaluated the suitability of adapting this logistic model to the 
actual GCC thickness/VFS in SAP distribution. 

Patients and Methods 

Subjects 

We retrospectively reviewed the medical records of patients with 
primary open-angle glaucoma who underwent both SAP by the 
Humphrey Visual Field Analyzer (Carl Zeiss Meditec, Dublin, 
CA) using the central 10-2 (HFA 10-2) program with a Swedish 
interactive threshold algorithm standard and macular GCC 
measurement by SD-OCT within a 3-month inter\'al at Keisho- 
kai lidabashi Eye Clinic, Tokyo, Japan. The HFA 10-2 inclusion 
period was from September 1, 2010 through August 31, 201 1. AU 
HFA 10-2 and macular GCC measurements by SD-OCT were 
performed during routine clinical practice. Glaucoma severity in 
the study population was determined by collecting the mean 
deviation (MD) values of HFA 24-2 within a year from the initial 
HFA 10-2 evaluation. Patients whose HFA 10-2 showed adequate 
reliability (fixation losses <33%, false positives and false negatives 
<20%) were included in the study. Analysis was performed only in 
the right eyes. All of the patients analyzed in the study underwent 
a complete ophthalmologic evaluation that included visual acuity 
testing, biomicroscopy of the anterior segment, intraocular 



pressure measurement by means of Goldmann applanation 
tonometry, and fundus examination. The inclusion criteria were 
a corrected visual acuity of 10/20 or better, a clear lens, and a 
normal retina. A 3-month testing window was selected to minimize 
the possibility of any glaucoma changes during the period between 
the structural and functional testing. 

The diagnosis of glaucoma was based on the evaluations of 
three glaucoma experts (DS, NO, and KY). Although observation 
of optic disc rim thinning, including an enlarged vertical cup to 
disc ratio (>0.7), or asymmetry (>0.2) between two eyes was used 
for diagnosing glaucoma, the experts also used their clinical 
expertise with or without these two criteria to make the diagnosis. 
Exclusion criteria included a history of other ocular or neurologic 
diseases that could affect the structural or functional measure- 
ments, intraocular surgery, and use of systemic or topical 
medications that might affect the perimetric results. 

Procedures followed the tenets of the Declaration of Helsinki, 
with the retrospecti\c protocol approved by the Institutional 
Review Board of Keisho-kai lidabashi Eye Clinic. Patient records 
and information were anonymized and de-identified prior to 
analysis. 

OCT measurements and GCC thickness analysis 

Macular GCC analysis was measured by SD-OCT (RS-3000; 
Nidek, Inc., Gamagori, Japan), which acquires 53,000 A-scans per 
second. The GCC map was centered on the fovea, with the central 
macula covered by square grid of 30°x30°. The GCC map was 
created by acquiring 128 vertical B-scans that were composed of 
512 A-scans. The 128 vertical B-scans were distributed in same 
interval within the 30° width. The GCC was automatically 
segmented from the internal Umiting membrane to the outer IPL 
boundary. 

The top panels of Figure 1 show representative data for the 
GCC thickness map analyzed by the RS-3000, the grayscale, and 
the SAP test point \ alucs d(;termined by the HFA 10-2 program. 
The built-in software divided the GCC map into three concentric 
circles (central, middle, and lateral circles), and placed horizontal 
and vertical lines across the fovea (Figure lA). The center of each 
circle was located at the fovea. The radii of the three concentric 
circles were 1.67°, 5°, and 10°, which respectively corrcspoiidccl to 
0.5, 1.5, and 3.0 mm in GuUstrand's schematic eye. The average 
thickness of the whole GCC scan field was measured as the mean 
GCC thickness of the area surrounded by the central and lateral 
circles. This division enabled the built-in software of the SD-OCT 
to measure the mean GCC thickness of the superior and inferior 
areas (Figure ID) or the eight segmented areas (Figure IE). The 
four areas surrounded by the c(;ntral and middle circles were 
defined as the 3-1 to 3-4 GCC-map divided area (GMDA), while 
the four areas surrounded by the middle and lateral circles were 
defined as the 6-1 to 6-4 GMDA (Figure IE). Images were 
excluded when the signal strength index was less than 7/10. When 
there were overt auto segmentation errors or obvious retinal 
abnormalities in the OCT, these images were also excluded. 

Adjustment of SAP test points to the GIVIDA 

To ensure the analysis between the GCC thickness and the VFS 
is correct, data adjustment is rc'(]uircd. Drasdo ct al. [12] rccendy 
reported finding relatively large displacement in six human retinas. 
Based on their ganglion cell receptive field model in the human 
visual field, they were able to calculate the predicted lateral 
displacement of the ganglion cells. This adjustment method has 
been described in detail elsewhere [7,11]. The predicted lateral 
displacement of the ganglion cells in the most centered visual field 
test points of the HFA 10-2 was 0.60 mm, while the displacement 
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Figure 1 . Methods of analysis for the macular ganglion cell complex thickness and standard automated perimetry. A) Representative 
ganglion cell complex (GCC) thickness map analyzed by RS-3000. Dividing lines and circles were drawn by the build-in software. B) Grayscale of the 
Humphrey visual field analyzer central 10-2 program in the same patient. C) Adjustment of the standard automated perimetry (SAP) test points to the 
GCC-map divided areas (GIVIDAs). Each inner GMDA (3-1 to 3-4) included three SAP test points (asterisks). Each asterisk of the SAP test points 
corresponds to the same color asterisk in the GCC map (top left). Each outer GMDA (6-1 to 6-4) includes 14 SAP test points. Note that the # is 
adjusted to the outer GMDA. D) Measurement of the GCC thickness in the superior and inferior areas. Numerical values in each area represent the 
mean GCC thickness (|im) of the area. E) Measurement of the GCC thickness in eight divided areas. These areas were named 3-1 to 6-4 GMDAs. 
Numerical values in each area represent the mean GCC thickness (|im) of the area. F) Mean visual field sensitivity is indicated in each GMDA. 
doi:1 0.1 371 /journal.pone.01 041 26.g001 



for the ganglion cells in tlie secondary centered visual field test 
points was 0.55 mm. The location of the most centered visual field 
test points and secondary centered visual field test points were 
0.48 mm and 0.94 mm, respectively. Thus, the eccentricity in 
GuUstrand's schematic eye for the ganglion cells corresponding to 
the most centered visual field test points was 1 .08 mm, while it was 
1.49 mm for the ganglion cells corresponding to secondary 
centered visual field test points. Based on these estimates, the 
inner area surrounded by the central and middle circles of the 
GCC map includes the four most centered visual field test points 
and the eight secondary centered visual field test points. Similarly, 
the outer area surrounded by middle and lateral circles of the 
GCC map includes another 56 test points (Figure lA, C). When 
the lateral displacement of the ganglion cells from the photore- 
ceptor cells in each quadrant is taken into consideration, the inner 
three SAP test points (asterisks in Figure IC) are included in the 
inner GMDA (i.e., 3-1, 3-2, 3-3, or 3-4). Similarly, adjustments of 
the outer 14 SAP test points, including the ^ point in Figure IC, 
led to their placement in the outer GMDA (i.e., 6-1, 6-2, 6-3, or 6- 



4). The mean VFS in each GMDA was calculated as the average 
of the adjusted SAP test point values (Figure IF). 

Multiple logistic model 

To investigate the relationship between the mean VFS and the 
explanatory variables, which included GCC thickness (xj: pan), age 
(X2: years), sex (xj: male = 0, female = 1), and spherical equivalent 
(SE) of refractive errors (x^: diopters), we adapted a multiple 
logistic model to the mean VFS and explanatory variable data in 
order to create formula 1 as follows: 



1 -f- exp(6o + bixi+ biXi + h^xi + b^xnO) q-j 
{y : mean VFS; d\, ho, hi, bj, 63, : unknown parameters) 



The above logistic regression analysis was solved using the 
Gauss-Newton method, with the estimators of the parameters 61, 
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and bo to b4 obtained byJMP 8 software (SAS, Gary, NC). Thep- 
values (po to ^4) for each explanatory variable parameter were also 
obtained. Adaptation between the real VFS in SAP and the 

estimated VFS from the multiple logistic model was analyzed using 
the test of goodness-of-fit [13]. This approaches use die idea of 
comparing the observed number of indi\ idu;ils with cacli outcome 
to the number expected based on the logistic regression equation 
[13]. These observed (0 = the real VFS in SAP) and expected 
(E = the estimated VFS from the logistic model) numbers are 
combined to form a goodness-of-fit chi square [13]. 

Large values of this test statistic indicate that the observed and 
expected results are very different, so that the regression equation 
has a poor fit to the data [13]. In contrast, small values of this test 
statistic indicate tliat the observed and expected r(;sults are similar, 
so the regression equation proves a good fit to the observations 
[13]. Alternatively, large P values (approaching 1) indicate a good 
fit and small P values (approaching 0) indicate a poor fit [13]. 

The distribution of the estimated VFS/GCC thickness dots 
from the model overlapped the real VFS/GCC thickness scatter 
pkjts. The significance of the model effect (i.e., predictability c)f the 
"VTS from explanatory variables) was examined by an analysis of 
variance (ANOVA), followed by the calculation of the coefficient 
of determination (83). 

Simple logistic model 

To construct a simple logistic model that depicts the approx- 
imate relationship between the GCG thickness and the mean VFS 
value, some of the other explanatory variables in formula (1) were 
removed in order to create formula (2), which is as follows: 

^ ei , . 

y l+exp[(xi-02)x(-(?3)] ^' 

The 02 parameter was defined as the GCG thickness at which 
the VFS decreased to half of the 9i, which represents the threshold 
of the obtained logistic curves. The 63 parameter reflects the 
steepness of the slope of the c:urves. Using the GCG thickness and 
the mean VFS data, the regression formula and parameters (Oj, 62, 
and 83) were obtained by the Gauss-Newton method. The 
adaptation between the real VFS in SAP and the estimated VFS 
from the simple logistic model was analyzed by the test of 
goodness-of-fit. The significance of the model effect was examined 
by an ANOVA, followed by the calculation of the R^. 
Subsequently, the logistic curve of the estimated VFS obtained 
from the simple logistic model/ GCG thickness was then drawn on 
the real VFS in the SAP/GGG thickness scatter plots. 

Statistical analysis 

Scatter plotting, Gauss-Newton method, test of goodness-of- 
fit, ANOVA, and other statistical analyses were performed using 
Excel 2007 (Microsoft, Redmond, WA) and JMP 8 software. A p 
value of <0.05 was considered statistically significant. 



Results 

Patient characteristics 

Out of the 95 candidate glaucoma eyes, 12 eyes were excluded 
due to overt segmentation error or other macular diseases such as 
idiopathic macular edema. As a result, the present study enrolled 
83 glaucomatous eyes (48 men and 35 women). The raw data of 
enrolled eyes are presented in Table SI. The mean age of the 
enrolled patients was 55.9±11.9 years, with a mean SE of 
-3.8±3.4 diopters, a MD of tiie HFA 24-2 of -6.5±7.2 dB and a 
MD of the HFA 10-2 of -5.99±6.4 dB. In 16 of the enrolled eyes, 
a glaucomatous visual field dc-fi-ct was not observed in the HFA 
24-2, even though the optic disc s in these eyes were diagnosed as 
glaucomatous. Figure 2 shows the histogram of the MD values of 
the HFA 24-2. There was no statistically significant relationship 
between any of the data pairs for age, SE and MD and there was 
no statistically significant difference for the age, SE, and MD 
between the men and women (data not shown). 

Multiple logistic model parameters 

Table 1 shows the 61 and the bo to b4 parameters, and the p- 
values for the explanatory variables obtained from the multiple 
logistic regression analysis. These values indicate the effect of the 
GCG thickness, age, sex, and SE of the refractive errors on the 
mean VFS in each area. The j&-values revealed that in all of the 
areas analyzed, the GCG thickness exhibited the greatest effect, 
and was the most statistically significant parameter among all of 
the selected explanatory variables. While the other explanatory 
variables were statistically significant in several areas, these effects 
were less than that observed for the GCG thickness. 

The estimators of h9, Iv? and b4 in the whole scan area, superior 
area, inferior area, and most of the GMDAs were all positive 
quantities, with some of their /^-values found to be statistit:ally 
significant. This suggests that patients who are older, have 
hyperopia, or are female might have less VFS as compared to 
patients who are younger, have myopia, or are male. While several 
estimators of b2, and b4 in GMDAs were inverted to negative 
quantities, none of their ^-values were significant. 

Fitting the multiple logistic model on the scatter plots of 
the VFS/GCC thickness 

Figure 3 shows the scatter plots between the real VFS/GCC 
thickness (blue circles). The distribution of the estimated VFS/ 
GCG thickness dots from the model (red squares) were found to 
overlap on the same scatter plots. 

Table 2 shows the /(-values from the test of goodness-of-fit, 
the F-and /(-values of the ANOVA, and the R^ that were obtained 
for each of the adapted multiple logistic model solutions. The x^ 
test of goodness-of-fit revealed that there was a good adaptation 
between the real VFS/GCC thickness and the estimated VFS/ 
GCG thickness from the multiple logistic model in all scatter plots 
except in the 6-3 areas. The F-and p-vakies from the ANOVA 
revealed that the multiple logistic models of all arenas w(-rc able to 
predict the VFS from the explanatory variables, with values 
statistically significant. 

Fitting the simple logistic model on the scatter plots of 
the VFS/GCC thickness 

Table 3 shows the 9i, 82, and 63 parameters obtained from a 
simple logistic model, which depicts the approximate relationship 
between the GCG thickness and the VFS. Table 4 shows the p- 
values of the x^ test of goodness-of-fit, the F- and p-values of the 
ANOVA, and the R^ that were obtained for the simple logistic 
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Figure 2. Histogram of glaucoma severity in the study population. The mean deviation values were measured by the central 24-2 program 
of the Humphrey visual filed analyzer. Patients with various stages of glaucoma were also included in the present study. 
doi:1 0.1 371/journal.pone.01 041 26.g002 



model solution. Logistic curves of the estimated VFS/GCC 
thickness from the simple logistic model (black curves) overlap 
the scatter plots of the real VFS/GCC thickness (Figure 3). 
Although the obtained logistic curves visually appeared to run 
through the center of the real VFS/ GCC thickness distribution on 
the scatter plots, the test of goodness-of-fit revealed that the 
estimated VFS from the simple logistic model faUed to adapt the 
VFS in SAP in several areas (3-3, 6-1, 6-3, 6-4). These results 
suggest that other explanatory variables, such as age, sex, and SE, 
were necessary for adapting the estimated VFS to the real VFS in 
SAP. On the other hand, the F- and j!>-values from the ANOVA of 
the simple logistic model revealed that the simple logistic models in 
all areas were able to predict the VFS from the GCC thickness. 
Even though these predictions were found to be statistically 



significant, overall these effects were less than those observed for 
the multiple logistic models (Table 4). 

Discussion 

Zeimer et al. [14] first reported that macular retinal thickness 
could be a criterion for the quantitative assessment of glaucoma. 
Axons, bodies, and dendrites of RGCs are believed to reside in the 
NFL, the GCL, and the IPL, respectively. Therefore, when 
detecting glauctjmatous changes, it may be more accurate to 
measure the thickness of these specific layers rather than the total 
macular thickness. Ishikawa et al. [3] developed a macular 
segmentation algorithm and used time-domain OCT to demon- 
strate its ability to quantify the glaucomatous change in the inner 
macular thickness. To detect glaucomatous changes in their study, 
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Figure 3. Visual filed sensitivity (VFS) versus ganglion cell complex (GCC) thickness. Scatter plots between the real VPS in standard 
automated perimetry and GCC thickness (blue circles) are shown. Distribution of the estimated VFS from the multiple logistic model/GCC thickness 
data overlapped on the same scatter plots (red squares). Curves of the simple logistic models (black curves) are also shown. 
doi:1 0.1 371 /journal.pone.01 041 26.g003 



they analyzed four retinal segments, the NFL, GCC, GCL+IPL 
and the total retina, and demonstrated that the macular GCC 
thickness provided optimal glaucoma detection of the highest 
repeatability. Since RS-3000 SD-OCT can also measure the 
thickness of GCL+IPL, we attempted to create scatter plots using 
the GCL+IPL thickness that corresponded to the VFS. However, 
since the data distribution was scattered too widely, we could not 
perform a proper analysis (data not shown). Thus, this is the 
reason why we selected the GCC thickness in our present analysis. 
Based on the previous findings, we speculate that auto-segmen- 
tation between NFL and GCL would be more difficult to perform 
than that between the ILM and vitreous cavity, even when using 
SD-OCT, including RS-3000. 

In the present study, the test of goodness-of-fit revealed that 
there was a good adaptation between the real VFS in SAP and the 
estimated VFS from the multiple logistic model in all of the scatter 
plots with the exception of the 6-3 area. In contrast, the estimated 
VFS from the simple logistic models were not able to adapt the 



VFS in SAP in a greater number of areas (3-3, 6-1, 6-3, 6-4). 
These results suggest that not only GCC thickness, but other 
explanatory variables, such as age, sex, and SE, are necessary for 
adapting the estimated VFS to the real VFS in SAP. Several 
studies have reported finding a structure and function relationship 
between the circumpapiUary retinal nerve fiber layer thickness and 
the perimetric sensitivity in glaucoma patients [15-19]. A number 
of studies have also examined the macular thicknesses of GCC or 
GCL+IPL, and the perimetric sensitivity [6-11]. Kim et al. [6] 
reported finding a relationship between the macular GCC 
thickness, the MD, and the visual field index of the HFA 24-2 
program. In their study, the relationships between the mean 
retinal nerve fiber layer/GCC thickness and the MD/visual field 
index were evaluated using regression analyses with a second- 
order polynomial formula. Furthermore, using a simple linear 
model based on previously published studies [7,16-18], Sato et al. 
[1 1] demonstrated there was a relationship between the relative 
field sensitivity and the GCL+IPL thickness in glaucoma. In line 
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Table 2. The outcome of the test of goodness-of-fit, ANOVA, and the of the multiple logistic models. 



Area 


test of goodness-of-fit 


ANOVA 








p-value 


A^value 


p-value 




whole field 


0.9966 


206.1 


<0.0001 


0.72 


superior area 


0.9979 


189.8 


<0.0001 


0.70 


inferior area 


0.3876 


267.3 


<0.0001 


0.77 


3-1 


0.9996 


179.1 


< 0.0001 


0.69 


3-2 


1 .0000 


180.5 


<0.0001 


0.69 


3-3 


0.1497 


408.4 


<0.0001 


0.83 


3-4 


0.9995 


252.6 


<0.0001 


0.76 


6-1 


0.2182 


160.4 


<0.0001 


0.66 


6-2 


0.9928 


290.7 


< 0.0001 


0.78 


6-3 


<0.0001 


179.4 


<0.0001 


0.69 


6-4 


0.7366 


236.2 


<0.0001 


0.74 
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with these previous reports, our current study also demonstrated 
that the effect of the GCC thickness on the VFS was greatest 
among all of selected explanatory variables. However, we also 
found that use of the single variable, GCC thickness, was 
insufficient and that other explanatory variables were necessary 
for properly adapting the estimated VFS to the real VFS. 

Our current study used a nonlinear multiple logistic model in 
order to avoid the limitations associated with the first-, second-, or 
higher-order polynomial models. These types of models cannot 
accurately present the so-called sigmoid or z-shaped curves that 
result from the plateau of the VFS or detectable limits of SAP. On 
the other hand, the multiple logistic model can portray such 
shapes, and was able to adapt well to our observed data for the 
test of goodness-of-fit. However, since the multiple logistic model 
was rejected in the 6-3 area, there could be further unknown 
explanatory variables that may increase the error for fitting models 
in the 6-3 area. Therefore, additional studies designed to look for 
other explanatory variables will need to be undertaken in the 
future. 

Previously there have been another report using models that 
GCC thickness reach a plateau when VFS were less than 20 dB 
[11]. This model stands on the assumption that there are layers 



that are not affected by glaucoma and the thickness of these layers 
are constant ( = GCC-thickness-plateau value) [20]. However, the 
GCC-thickness-plateau model has a severe defect that the curve of 
the model does not cover the GCC-thickness range less than 
GCC-thickness-plateau value. So in the case that GCC-thickness 
was less than GCC-thickness-plateau value, expected VFS value in 
this model would be mandatory regarded as 0 db, which might 
cause severe errors between expected and observed VFS values. 
Same problem would occur whenever certain lower limit of GCC- 
thickness values was placed in any models. So we avoid placing 
any limitation in GCC-thickness range so that expected VFS 
values would be obtained close to observed VFS values in any 
GCC-thickness values. The gradual slope of our simple and 
multiple logistic models in lower GCC-thickness range is the 
original merit, since it works to represent expected values of VFS 
at any GCC-thickness values. When compared presented simple 
logistic model and previous models, values of presented simple 
logistic model exceeded that of previous models (Table 5), 
suggesting that presented logistic models functions have better 
predictability of the VFS from explanatory variables compared to 
linear model, second-order polynominal, third-order polynominal 
or GCC-thickness-plateau model. 



Table 3. The estimated parameters of the explanatory variables obtained from the 


simple logistic regression analysis. 












whole field 


33.479 


56.629 


0.080 


superior area 


34.303 


54.999 


0.076 


inferior area 


31.781 


58.076 


0.139 


3-1 


34.239 


46.506 


0.082 


3-2 


35.114 


51.061 


0.057 


3-3 


33.386 


59.917 


0.135 


3-4 


32.656 


59.154 


0.197 


6-1 


33.758 


49.706 


0.099 


6-2 


32.765 


58.038 


0.088 


6-3 


32.129 


52.983 


0.212 


6-4 


31.490 


56.007 


0.102 
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Table 4. The outcome of the test of goodness-of-fit, ANOVA, and the of the simple logistic models. 



Area 


test of goodness-of-fit 


ANOVA 








p-value 


A^value 


p-value 




whole field 


0.9606 


144.8 


<0.0001 


0.64 


superior area 


0.9344 


128.3 


<0.0001 


0.61 


inferior area 


0.0698 


210.1 


<0.0001 


0.72 


3-1 


0.9990 


115.7 


< 0.0001 


0.59 


3-2 


0.9501 


91.0 


<0.0001 


0.53 


3-3 


0.0413 


295.7 


<0.0001 


0.78 


3-4 


0.9982 


234.4 


<0.0001 


0.74 


6-1 


0.0317 


110.1 


<0.0001 


0.58 


6-2 


0.0593 


104.2 


< 0.0001 


0.56 


6-3 


<0.0001 


156.3 


<0.0001 


0.66 


6-4 


0.0084 


106.2 


<0.0001 


0.57 
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In general, the logistic regression model is used for the 
dichotomous distribution. Concerning propriety of using logistic 
models for VFS/GCC thickness distributions, we consider VFS 
reflects the survival rate of functioning ganglion cells, which is 
probability of dichotomous variation, since VFS consists of sum of 
responses of the functioning ganglion cells. When total ganglion 
cells are healthy and functioning properly, VFS would be 
maximum value ( = 0i). As more and more ganglion cells suffer 
or die from glaucoma and lost their functions, survival rate of 
functioning ganghon cells decreases, and VFS/Oi rate also 
decreases at the same time. When the assumption stands that 
VFS consists of sum of responses of the fvmctioning ganglion cells, 
the survival rate of functioning ganglion cells and VFS/9i rate 
would be correspond, and VFS/Qi would also be regarded as 
probability of dichotomous variation. We speculate this is why 
logistic models can be adapted to VFS/GCC thickness distribu- 
tions so good. Unfortunately, it is just a speculation since there are 
no ways to measure the number of functioning ganglion cells in 
vivo. Further evaluations about relationship between VFS and 
survival rate of functioning ganglion cells are needed, and a 
method for measuring functioning ganglion cells is to be improved. 

As the GCC thickness of the nasal area includes certain portions 
of the optic nerve fibers that originate from the ganglion cells in 
the temporal area, we were initially afraid that the fitness of the 



simple or multiple logistic models to the nasal area might be 
statistically rejected. However, our results showed these models 
were statistically accepted in the nasal areas. We speculate that 
these results are related to the anatomical characteristics of the 
retinal nerve fibers. Since the normal NFL thickness in the macula 
has been reported to be about one-half of the GCL-HiPL thickness, 
the influence of the NFL thickness may not be great enough to 
increase the error of the fitness [4] . 

Since our multiple logistic model included a greater number of 
explanatory variables, its predictability was more accurate than 
that of the simple logistic model. However, the simple logistic 
model had some merits of its own, as it can visualize the logistic 
sigmoid curve that depicts the rough relationship between the 
GCC thickness and the mean VFS value. In addition, the 
thresholds of the obtained logistic curves were simply represented 
as the parameter Gj, which may be much easier to image than the 
multiple logistic model. 

The analysis by these models demonstrated several difierent 
findings. First, multiple regression analysis revealed that patients 
who are older, have hyperopia or are female might have less VFS. 
Most previous studies that have examined the structure and visual 
field sensitivity relationship in glaucoma did not include such 
explanatory variables. Similarly, as our test of goodness-of-fit 
was rejected in the 6-3 area, our current study also may not have 



Table 5. Comparison of values between presented models and previous models. 





whole field 


superior area 


inferior area 


Simple logistic model (present study) 


0.64 


0.61 


0.72 


Multiple logistic model (present study) 


0.72 


0.70 


0.77 


Linear (Kim Na, et al [6])* 


0.26 


n/a 


n/a 


Second-order Polynominal (Kim Na, et al [6]}* 


0.30 


n/a 


n/a 


Third-order Polynominal (Kim Na, et al [6])* 


0.31 


n/a 


n/a 


GCA-thickness-plateau model (Sato S, et al [11])^ 


0.38 


0.24 


0.45 



GCA = ganglion cell layer and inner plexiform layer. 

* values between ganglion cell complex thickness and mean deviation value of 24-2. 
^R^ values between GCA thickness and mean visual filed sensitivity by microperimetry. 
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fully included all of die explanatory variables. Therefore, die 
mechanisms responsible for die effect of these explanatory 
variables on the VFS remain to be elucidated. 

Since the plateau of the logistic curves suggests that the GCC 
thickness change in glaucoma occurs much earlier than the SAP 
sensitivity loss, SD-OCT measurements of the GCC thickness can 
be used to detect early glaucomatous damage. Thus, GCC 
thickness may be one of the most advantageous measurements for 
preperimetric or early perimetric glaucoma, which can be dilRcult 
to evaluate when using SAP. 

In conclusion, after assessing explanatory variables such as the 
GCC thickness and examining the SAP sensitivity in patients with 
glaucoma, we were able to create a logistic model that can be used 
to depict the relationship between these variables. With the 
exception of the nasal-inferior-outer area (6-3), all of the multiple 
logistic models exhibited a good ability to fit the real GCC 
thickness/VFS in the SAP distribution to the test of goodness- 
of-fit. Our ANOVA results also showed that the models exhibited 
a statistically significant ability to predict the VFS from 
explanatory variables, which suggests that the model is a reliable 
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